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Integrating  Numerical  Models  with  Data  Analysis  in  Site  Assessment 

Timothy  R.  Keen  (keen@nrlssc.navy.mil)  and  James  D.  Dykes 
(Naval  Research  Laboratory,  Stennis  Space  Center,  Mississippi) 

ABSTRACT:  The  characterization  and  assessm  ent  of  contaminated  sites  can  be  m  ade 
more  effective  with  the  interactive  use  of  in  situ  data,  historical  records,  and  results  from 
predictive  process-based  numerical  models.  These  three  components  can  be  difficult  to 
merge  because  of  space  and  tim  e  scales,  as  well  as  data  form  ats  and  availability.  We  are 
developing  an  Application  and  Programming  Interface  (API)  tool  to  make  this  procedure 
more  accessible  to  a  ran  ge  of  interested  groups  .  This  tool  consists  of  two  com  ponents: 
(1)  a  GUI-based  toolbox  of  data  processing  and  analysis  software  and  (2)  a  num  erical 
modeling  system  for  hydrodynam  ics,  sedimentation,  and  m  ass  transport.  The  toolkit 
includes  modules  to  analyze  data  in  four  dim  ensions  and  merge  different  data  types  as 
well  as  compute  deriv  ed  variables.  The  m  odeling  system  includes  models  that  calculate 
waves;  circulation  due  to  tides  and  the  wi  nd;  entrainment  of  bottom  sediment;  and  the 
transport  of  dissolved  and  suspended  m  aterial.  The  com  ponents  of  the  API  have  been 
tested  on  the  potential  entrainm  ent  and  transport  of  contaminants  in  San  Francisco  Bay, 
California,  and  St.  Louis  Bay,  Mississippi.  The  models  have  proven  valuable  at  predict¬ 
ing  the  movement  of  dissolved  and  suspended  contaminants  both  within  enclosed  bays  as 
well  as  in  estuaries  and  the  inner  shelf.  The  results  indicate  that  integration  of  the  models 
with  available  data  for  planning  containment  and  remediation  strategies  would  reduce  the 
effort  required  to  characterize  a  contaminated  site  using  field  methods  alone. 

INTRODUCTION 

The  U.S.  Navy  has  m ore  than  40  coastal  and  nearshore  contaminated  sediment  sites 
that  are  included  on  the  National  P  riorities  List  (Superfund  sites)  (CNO,  2005).  Navy 
remediation  policy  requires  that  the  source(s)  for  these  sites  must  be  identified  and  that  a 
monitoring  plan  must  be  in  place  before  data  collection  begins.  Both  of  these  tasks  can  be 
made  easier  by  the  use  o  f  mathematical  models  of  the  site-specific  processes  that  impact 
contaminant  behavior  and  public  health  risks.  For  example,  inactive  contaminant  sources 
can  be  identified  from  historical  records  and  the  current  distribut  ion  of  sediment-bound 
contaminants,  but  this  straightforward  approach  cannot  predict  remobilization  by  physical 
mechanisms. 

Contaminant  transport  by  water  and  sedim  ents  can  be  simulated  in  numerical  hydro- 
dynamic,  sediment  dynamics,  and  water  quality  models.  These  m  odels  are  capable  of 
simulating  post-depositional  processes  and  can  help  identify  sources  and  evaluate  remo¬ 
bilization.  The  quantitative  pred  ictions  from  numerical  models  can  be  used  to  define 
zones  of  potential  future  contam  ination  and  assist  in  determining  an  efficient  sam  pling 
strategy  for  monitoring.  Sediment  and  contaminant  transport  modeling  is  often  conducted 
on  a  site-specific  basis.  This  approach  works  especially  well  for  rivers  and  lakes  (Con¬ 
nolly  et  ah,  2000)  but  accuracy  dep  ends  on  dense  sampling  because  local  models  cannot 
include  the  impacts  of  external  influences.  With  the  recognition  of  pollution  problems  in 
estuaries,  the  simulation  of  contaminant  dispersal  in  open  bodies  of  wa  ter  has  become 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
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more  widespread.  The  accurate  prediction  of  contaminant  transport  in  coastal  waters 
requires  comprehensive  hydrodynamic  descriptions  (i.e.,  river  input,  tides,  wind-driven 
flow,  heat  and  m  ass  mixing,  and  surface  waves)  .  The  U.S.  EPA  encourages  the  use  of 
numerical  models  for  pollutant  transpor  t  (e.g.,  the  Chesapeake  Bay  Program  :  www. 
chesapeakebay.net),  which  is  especially  im  portant  for  long-term  effects  caused  by  rare 
events  (e.g.,  severe  storms)  and  the  cumulative  effects  of  weaker  processes. 

Numerical  models  are  typically  used  only  where  there  is  public  interest  or  a  potential 
health  threat  (e.g.,  Aldridge  et  al.,  2003).  NAVFAC’s  Implementation  Guide  recommends 
that  modeling  should  be  considered  as  a  tool  to  integrate  site  data  (NAVFAC,  2003),  but 
conceptual  site  m  odels  are  routinely  used  at  Navy  sites  instead  of  numerical  models. 
Nevertheless,  the  Naval  Oceanographic  Office  (NAVOCEANO)  and  Naval  Research 
Faboratory  (NRF)  have  developed  and  successf  ully  validated  several  models  for  littoral 
environments  (Harding  et  al.,  1999). 

INTEGRATED  AQUATIC  DATA  ANALYSIS  AND  SIMULATION 

Our  approach  consists  of  two  com  ponents:  (1)  a  suite  of  num  erical  models  and  the 
software  required  for  their  use  and  (2)  an  analysis  toolkit  built  on  Geographic  Informa¬ 
tion  System  (GIS)  lib  raries.  This  system  is  being  actively  developed  and  the  results 
presented  herein  represent  an  in  termediate  level  of  implementation  rather  than  a  com¬ 
pleted  API. 

Numerical  Models.  Numerical  models  are  routinely  used  for  atm  ospheric  and  oceano¬ 
graphic  forecasting  at  production  centers,  forw  ard  centers,  and  even  at  the  fleet  level 
(Flather,  2000;  Burnett  et  al.,  2001).  The  fo  rward  deployment  of  ocean  modeling  is  sim¬ 
plified  somewhat  from  the  produ  ction  centers  becaus  e  of  reduced  com  putational 
resources,  but  local  models  are  su  pplied  with  required  boundary  conditions  from  the 
major  centers.  This  approach  has  been  dem  onstrated  for  different  physical  processes  in 
several  regions  (Keen  and  Holland,  2010). 

Simulations  of  mass  transport  in  e  stuaries  and  nearshore  areas  are  d  ependent  on 
water  levels  and  currents  calculated  by  circulation  models.  Computing  chemical  behavior 
is  further  reliant  on  salinity  and  temperature.  Focal  hydrodynamic  models  are  more  flexi¬ 
ble  than  global  and  basin  m  odels  because  they  do  not  incorporate  d  ata  assimilation 
schemes.  A  tailored  model  can  thus  be  implemented  for  a  specific  problem.  In  addition  to 
flow  calculations,  it  is  o  ften  necessary  to  simulate  waves  and  mass  transport  processes 
(e.g.,  sedimentation),  requiring  either  additional  num  erical  models  or  com  prehensive 
observations.  The  examination  of  complex  sites  where  multiple  sources  are  present  and 
remobilization  is  possible  would  likely  requ  ire  the  simulation  of  the  wind,  waves,  cur¬ 
rents,  temperature  and  salinity,  sedim  entation,  boundary  layer  processes,  and  m  ass 
transport  within  a  sm  all  area,  as  well  as  the  surrounding  estuary  or  ocean.  Re  mediation 
experts  thus  find  them  selves  forced  to  choos  e  between  using  a  param  etric  model  that  is 
dependent  on  random  sampling,  and  using  num  erical  models  that  require  expert  knowl¬ 
edge  for  reliable  operation. 

One  solution  to  this  problem  is  to  implement  numerical  models  of  physical  processes 
in  the  littoral  ocean  in  the  following  manner.  First,  it  is  a  simple  matter  to  take  advantage 
of  operational  ocean  models  run  by  NAVOCEANO.  These  fields  of  water  levels,  currents 
and  waves  can  be  a  first  approxim  ation  of  environmental  forcing  in  m  any  areas.  The 
second  step  is  to  utilize  local  models  in  a  reach-back  mode,  which  is  largely  implemented 


to  provide  environmental  support  for  Navy  operations  (Malley  et  al.,  2005).  Finally,  this 
approach  can  be  im  plemented  as  a  network-  centric  system  as  has  been  proposed  for 
operational  forecasting  (Jimenez,  2009). 

Analysis  Toolkit.  NAVOCEANO  has  im  plemented  an  ope  rational  support  center  for 
antisubmarine  warfare  (ASW)  oceanographic  modeling  and  forecasting.  This  reach  -back 
cell  (RBC)  is  staffed  by  subject-m  atter  experts  whose  analyses  are  aided  by  the  AS  W 
Reach-Back  Cell  Oceanography  Analysis  System  (ARCOAS).  ARCOAS  is  a  GIS-based 
toolkit  incorporating  a  custom  ized  ArcMAP  application  tailored  for  meteorological  and 
oceanographic  (METOC)  analysis  of  data  a  nd  model  forecasts.  ARCOAS  has  extensive 
capabilities,  including  georeferencing  using  a  standard  datum  (e.g.,  NAD83),  the  statisti¬ 
cal  analysis  of  different  data  bases  on  their  o  riginal  locations  or  co  Hated  for  sp  ecified 
times,  spatial  and  temporal  averaging,  and  the  inclusion  of  other  data  (e.g.,  USA  Prime 
Imagery  server). 

ARCOAS  can  easily  be  expanded  to  in  elude  additional  M  ETOC  and  other  geo¬ 
physical  tools  as  the  need  aris  es  because  it  is  programm  ed  using  an  inherently  flexible 
and  full-featured  platform.  This  advantage  pe  nnits  access  and  analys  is  of  all  data  and 
information  in  one  place.  ARCOAS  both  addresses  operational  Navy  requirem  ents  and 
lends  itself  to  interoperability  with  other  systems  using  similar  protocols  and  formats.  A 
primary  objective  is  to  provide  critical  enviromnental  characterization  of  the  battle  space 
within  the  c  oncept  of  network-centric  operations  using  a  se  rvices-oriented  architecture 
(SOA)  (Meyer,  2007).  In  line  with  this  Navy  METOC  requirement,  the  RBC  co  ncept 
provides  the  ability  for  forward-deployed  team  s  to  reques  t  in-depth  analyses  from  the 
RBC,  thus  freeing  them  to  provide  tailored  support  to  th  eir  customers.  Some  ARCOAS 
tools  may  also  be  suited  for  users  outside  of  the  RBC. 

EXAMPLE  PROBLEM:  HUNTERS  POINT  SHIPYARD 

The  south  basin  at  Hunters  Point  Sh  ipyard  (SBHP)  is  a  small  (<  0.5  km2)  tidal  basin 
in  San  Francisco  Bay  (SFB)  (Figure  1).  Pote  ntial  remobilization  in  SBHP  depends  on 
whether  sediment-bound  contaminants  can  escape  into  SFB  by  natural  m  echanisms.  To 
address  this  problem ,  it  is  necessary  to  exam  ine  the  potential  release  o  f  contaminants 
from  the  bottom  and  to  evaluate  their  poten  tial  transport  by  physical  processes.  This 
complex  problem  has  been  studied  using  wind  forecasts  from  the  Navy  Operational 
Global  Atmospheric  Prediction  System  (NOGAPS),  and  Pacific  Ocean  tern  perature, 
salinity,  and  sea  surface  height  (SSH)  from  the  global  Navy  Coastal  Ocean  Model 
(NCOM).  These  products  supply  boundary  conditions  to  local  grids  for  NCOM  to  calcu¬ 
late  currents  and  SSH,  and  SWAN  (Si  mulating  Waves  Nearshore),  which  computes 
significant  wave  height  (SW  H).  The  predicted  curren  ts  and  waves  are  used  to  com  pute 
sediment  entrainment  with  the  Litto  ral  Sedimentation  Model  (LSM).  The  currents  are 
also  used  to  calculate  dissolved  tracer  concentrations  and  particle  transport. 

Modeled  Hydrodynamics  at  Hunters  Point.  This  example  focuses  on  the  contam  inant 
remobilization  during  periods  of  southeasterl  y  winds  and  tides.  Th  e  largest  waves  in 
SBHP  occur  when  the  wind  is  from  the  SE  due  to  the  increased  fetch  (Ziimnerman  et  al., 
2008;  Keen  and  Holland,  2010).  The  NOGAPS  wind  is  southeasterly  (m ax  =  13  m/s) 
between  6  and  10  January  2004.  The  predicte  d  SWH  from  SWAN  in  SBHP  exceed  s 
35  cm  on  8  January.  The  ebb-tide  model  curre  nts  (short-dash  line  in  Figure  2)  reach 


6  cm/s  toward  the  SE  along  the  eastern  shoreline  of  SBHP  despite  the  southeasterly  wind. 
This  flow  extends  from  the  NE  comer  of  the  inner  basin  to  the  outer  basin.  The  tidal  flow 
is  into  SBHP  during  the  flood  tid  e.  This  inflow/outflow  cycle  is  repeated  daily.  There  is 
also  a  strong  ebb  tide  flow  (>  20  cm/s)  from  Yosemite  Creek  (long-dash  line  in  Figure  2) 
in  the  NW  end  of  SBHP.  The  model  currents  during  high  water  show  a  counterclockwise 
transport  within  SBHP  that  is  in  agreem  ent  with  the  observations  and  residual  sediment 
fluxes  (Zimmerman  et  al.,  2008).  This  transport  regime  includes  a  SE  flow  driven  by  the 
ebb-tide  that  the  observations  from  2001  did  not  capture.  The  model  also  reveals  a  wide¬ 
spread  flow  into  the  outer  basin  and  transport  along  the  eastern  shore. 
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FIGURE  1.  Grids  used  for  numerical  simulations  in  this  study: 

(A)  the  SFB  grid;  (B)  the  Hunters  Point  (HP)  grid;  and  (C)  the  SBHP  grid. 

The  box  in  A  indicates  the  domain  of  B,  and  the  box  in  B  is  the  SBHP  grid. 

Resuspension  of  Bottom  Sediment.  The  predicted  max  shear  stress  rm  is  8.7  Pa.  The 
depth  of  erosion  can  thus  be  qualitatively  compared  to  the  resu  Its  from  laboratory  ex¬ 
periments  on  sediments  from  Hunters  Point  (Zimmerman  et  al.,  2008).  The  Sedfl  ume 
device  was  used  to  m  easure  the  erosion  rates  for  homogenized  sediments  after  1  hr  and 
16  d  of  consolidation.  The  erosion  rate  for  se  diments  at  1  and  13  cm  in  a  core  (16  d 
settling)  at  a  shear  stress  of  6.4  Pa  was  0.01  and  0.008  cm/s,  respectively.  The  same  core 
had  a  constant  erosion  rate  of  approxim  ately  0.03  cm/s  at  a  shear  stress  of  8  Pa.  This 
result  suggests  that  the  surface  sediments  would  have  been  instantly  removed. 


FIGURE  2.  Summary  of  currents  and  tracer  predictions  for  SE  wind  event  during  a  tidal 
cycle  on  8-9  January.  The  maximum  tidal  outflow  is  20  cm/s  (long  dashed  line)  and  tidal 
inflow/outflow  is  6  cm/s  (short  dashed  line).  The  wind  speed  is  labeled  and  the  direction  is 
indicated  by  the  arrow.  The  tracer  concentration  is  contoured  after  24  h  of  release  from  the 
locations  indicated  by  stars  (dimensionless).  The  measured  concentration  of  PCBs  (pg/kg) 
is  indicated  by  the  solid  circles  (Battelle  et  al.,  2002). 

The  LSM  calculates  total  suspended  solids  (TSS)  as  well  as  an  activ  e  layer  depth, 
which  is  the  thickness  of  disturbed  sediment,  and  the  equiva  lent  bed  depth  of  sedim  ent 
particles  in  suspension.  The  predicted  TSS  (Figure  3)  reaches  peaks  of  450  and  360  mg/L 
during  consecutive  low  tides;  the  resulting  active  layer  exceeds  14  cm  and  the  resuspen¬ 
sion  depth  is  0.54  mm.  The  sim  ilarity  of  the  measured  (400  to  600  mg/L  in  2001)  and 
predicted  TSS  suggests  an  erosion  depth  of  a  bout  10  cm,  below  whi  ch  the  sedim  ent 
would  have  been  too  consolidated  to  erode. 

The  depth  of  erosion  can  be  used  to  evaluate  the  release  of  contam  inants  from  the 
bottom  sediment.  For  example,  the  measured  concentrations  of  metals  from  grab  samples 
(Table  B3  in  Battelle  et  al.,  2002)  suggest  el  evated  levels  within  the  predicted  erosion 
depth.  These  metals  could  thus  be  introduced  to  the  water  column  during  SE  wind  events. 

The  mechanism  for  their  exchang  e  between  the  sediment  and  water  would  need  to  be 
represented  by  a  water  quality  m  odel  that  is  not  part  of  our  system .  The  available  profile 
data  for  SBFIP  are  not  easily  analy  zed  because  typical  sample  intervals  were  2  feet  (60 
cm)  except  for  Pb210  data,  which  were  collected  at  ~10  cm  intervals. 


FIGURE  3.  Hourly  TSS  20  cm  above  the  bed  computed 
by  LSM  in  SBHP  (water  depth  =  1.82  m). 


Contaminant  Transport.  Contaminants  within  the  seafloor  can  be  released  when  the 
bottom  sediment  is  disturbed.  The  timing  of  the  potential  releases  in  this  study  is  directly 
related  to  the  wave  energy  predicted  by  SWAN  through  the  bottom  shear  stress.  The  most 
likely  times  for  contaminants  to  be  released  into  the  tidal  basin  are  when  TSS  is  high 
(Figure  3).  This  potential  release  is  sim  ulated  in  the  m  odel  using  a  concentration  of  5 
(units  are  arbitrary)  at  the  locations  indicated  in  Figure  2. 

The  initial  tracer  release  coincides  with  the  first  resuspension  event,  which  occurs  at 
low  tide.  During  the  ensuing  flood-tide,  the  tracer  released  within  the  inner  basin  is  trans¬ 
ported  westward  and  collects  along  the  wester  n  shore.  The  CCW  circulation  within  the 
south  basin  causes  accumulation  of  the  tracer  al  ong  the  western  side  of  the  inlet  as  well. 
The  flood  tide  works  with  the  SE  wind  to  tran  sport  the  surface  water  back  to  the  NW  and 
into  the  inner  basin.  The  following  ebb  tide  transports  the  tracer  along  the  eastern  shore 
into  the  outer  basin  (Figure  2). 

A  second  approach  to  evaluating  the  paths  of  contaminants  released  from  sediments 
at  SBHP  is  to  use  particles  as  Lagrangian  tracers.  This  m  ethod  is  useful  for  complex 
flows  where  the  kinem  atics  of  particle  beha  vior  are  not  well  constrained.  The  particle 
model  uses  the  currents  calculated  by  NCOM  on  the  SBFIP  grid  for  January  2004  (Figure  4). 
After  31  da  ys  of  transport,  the  h  eavier  particles,  which  ar  e  representative  of  silt,  are 
primarily  distributed  around  Hunters  Point  (yellow  points  in  Figure  4)  and  becom  e 
trapped  within  the  southern  basin.  The  lighter  particles  (aqua  circles)  have  been  trans¬ 
ported  into  SFB.  The  interes  ting  result  from  this  simulation  is  that  particles  with  finite 
settling  velocities  have  reached  the  m  ain  bay.  This  tran  sport  occurs  over  several  days 


under  a  range  of  flow  conditions.  Th  e  implication  is  that,  even  if  toxins  adsorbed  to  clay 
particles  require  several  days’  exposure  to  sea  water  to  become  dissolved,  it  is  likely  that 
some  contaminant  would  be  dissolved  into  the  waters  of  the  main  bay. 


FIGURE  4.  Particle  distributions  near  Hunters  Point  computed  by  a  particle  tracking 
algorithm  for  January  2004.  Particles  were  released  at  the  locations  indicated  by  the 
circles.  The  aqua  particles  have  a  settling  speed  of  10'  mm/s,  which  is  representative 
of  clay  floes.  The  yellow  particles  have  a  settling  speed  of  10'2  mm/s,  which  is  more 
typical  for  large  floes  and  silt.  The  bathymetry  is  represented  by  shading.  The  vectors 
are  the  surface  currents  calculated  by  NCOM  on  the  30-m  grid.  Current  speed  (m/s)  is 
indicated  by  shading  as  shown  on  the  colorbar.  The  inset  map  shows  the  particles  after 
15  days  (red  pluses)  and  28  days  (blue  pluses)  for  February  2004.  The  currents 
are  from  the  NCOM  SFB  grid.  One  particle  was  released  every  24  hours. 


SUMMARY 

The  paper  has  two  objectives:  (1)  shed  light  on  the  question  of  whether  and  how 
sediment-bound  toxins  in  a  sm  all  tidal  basin  can  escape  into  San  Francisco  Bay  and  (2) 
briefly  demonstrate  the  potenti  al  integration  of  numerical  modeling  methods  with  data 
collected  from  the  field  in  studying  this  problem .  Part  (1)  is  a  straightforward  application 
of  multiple  models  to  the  prob  lem  but  it  is  complicated  by  the  importance  of  environ¬ 
mental  factors  that  op  erate  at  a  range  of  scales.  This  problem  demonstrates  the 
importance  of  tides  and  atmospheric  events  like  the  southeasterly  winds  in  SFB.  These 
environmental  factors,  which  are  external  to  the  study  area,  are  critical  in  understanding 
the  hydrodynamics  and  potential  release  of  sediment-bound  contaminants.  The  internal 
factors  like  local  wind  and  still-water  depth  were  not  as  criti  cal  in  this  example  but  they 
play  important  roles  for  the  long-term  stabili  ty  of  the  bottom  sedi  ment.  There  is  no  easy 
way  to  include  both  internal  and  external  environmental  factors  in  a  numerical  study  but 
this  example  shows  one  approach  that  does  not  require  the  development  of  new  models. 
It  does,  however,  rely  on  access  to  operational  or  archived  environmental  forcing.  This  is 
where  the  reach-back  cell  concept  becomes  important,  as  has  been  the  case  for  ASW. 

This  study  has  dem  onstrated  some  of  the  problems  associated  with  multi-s  cale, 
multi-physics  simulations.  As  with  most  estuary  models,  the  water  depth  is  critical.  The 
shallow  water  in  the  so  uth  basin  is  a  dete  nnining  factor  in  sedim  ent  resuspension.  A 
small  uncertainty  in  depth  can  decrease  the  wave-current  bottom  stress  below  the  critical 
threshold  for  entrainment  and  complete  ly  change  the  r  esults.  Thus,  it  is  im  portant  to 
simulate  the  tidal  amplitude  and  phase  accurately  so  that  the  water  depth  is  correct.  The 
wind  direction  was  critical  in  this  exam  pie  because  of  the  generation  of  larger  waves 
from  the  SE.  The  range  of  potential  environm  ental  factors  that  can  im  pact  the  hydrody¬ 
namics  and  sediment  dynamics  in  an  estuary  is  a  strong  m  otivation  for  implementing  the 
RBC  concept  using  network-centric  software  like  ARCOAS.  This  approach  can  bring  the 
computational  and  data  resources  of  the  enti  re  defense  establishment  into  the  hands  of 
remediation  managers  at  forward  deployments. 
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